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Q_) We show that extraordinary magnetoresistance (EMR) arises in systems consisting of 

C^ two components; a semiconducting ring with a metallic inclusion embedded. The im- 

,__! portant aspect of this discovery is that the system must have a quasi-two-dimensional 

character. Usi ng the same materials and geometries for the samples as in experiments by 

SSolin et o/LLLil we show that such systems indeed exhibit a huge magnetoresistance. The 
magnetoresistance arises due to the switching of electrical current paths passing through 
^-J the metallic inclusion. Diagrams illustrating the flow of the current density within the 

samples are utilised in discussion of the mechanism responsible for the magnetoresistance 
effect. Extensions are then suggested which may be applicable to the silver chalcogenidcs. 
I Our theory offers an excellent description and explanation of experiments where a huge 

' {J) magnetoresistance has been discoveredHEl 

a 

Q Keywords: Magnetoresistance; Current paths; Strongly inhomogencous. 

1. Introduction 

The magnetoresistance effect known as extraordinary magnetoresistance (EMR) was 

coined by Solin et alW after the experimental discovery in 2000. This effect showed 

_i_ extremely large room temperature magnetoresistance values of a million percent 

l/^ in a 5T applied magnetic field-M These values were found in composite van der 

r — Pauw disks consisting of a conducting and semiconducting region of non-magnetic 

materials. As a result of the observed properties of EMR it has been proposed that 

improvements could be achieved in the read heads of magnetic disk driveSH 

In 1997 Xu et al. 3 discovered a large magnetoresistance in non-magnetic silver 
chalcogenides Ag 2+s Se and Ag 2+s Te. The addition of a small excess amount (5) of 
silver atoms into the structure of the semiconductors Ag 2 Se and Ag 2 Te (with no 



X 

appreciable magnetoresistance themselves) caused a large linear magnetoresistance 
to be observed. At room temperature and in a magnetic field of 5.5T the magnetore- 
sistance was 200% with no sign of saturation. The magnetoresistance displayed a 
linear dependance on magnetic field which remained in low fields^. A model for this 
magnetoresistance has been proposed 6 ! 7 ! 8 ! which considers the silver chalcogenides 
as consisting of two components: a semiconductor and a conducting component. 
This is analogous to the EMR effect in systems where it may arise. Bulgadaev 
and Kusmartse\12112l have found explicit expressions for the magnetoresistance of 
strongly inhomogeneous two-phase systems by developing a method of conformal 
mapping transformations. This method utilises the exact dual transformation in re- 
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lating the effective conductivity of plannar inhomogeneous two-phase systems with 
and without an applied magnetic field. They have studied three models: the random 
droplet model (RDM); the random parquet model (RPM); and an effective medium 
model (EMM). 

Parish and Littlewood 11 have modelled a strongly inhomogeneous conductor us- 
ing a random resistor network made from four terminal resistors, and have shown 
that the properties of the magnetoresistance are similar to those of the silver chalco- 
genides. Abrikosov 12 proposed an explanation of the abnormal magnetoresistance 
observed in the silver chalcogenides, namely quantum magnetoresistance, based on 
the assumption that in such systems gapless dirac fermions arise. Recent experimen- 
tal work has verified that the magnetoresistance effect in the silver chalcogenides 
is somehow related to the microstructure of the material, which has been identified 
and investigated 13 l 14 l 

In this paper van der Pauw disks consisting of a conducting and semiconducting 
region of non-magnetic materials are considered. Their electrical properties, such 
as conductivity and the distribution of the electric field and current, are studied 
using finite clement simulations. These were based on the EMR geometry used 
by Solin et alW when the effect was first discovered. The experimental results are 
verified and the mechanism for the effect discussed with the use of streamline plots 
of total current density. The paper is structured as follows: firstly, a description of 
the approach taken in order to produce the simulations including deviations from 
the experimental geometry; secondly, the magnetoresistance results are given along 
with the streamline plots; and finally, discussion of the results in comparison with 
previous experiments, with the mechanism for the magnetoresistance offered as well 
as proposals for extension. 

2. Simulation approach 

We modelled the experimental construction of the composite van der Pauw disk, 
as used by Solin et aZP, in the discovery of the EMR effect. The experimental 
construction consisted of a circular disk of Indium Antimonide (InSb) of radius r& = 
0.5mm with a concentric inclusion of Gold (Au) of radius r a . Different geometries 
were achieved by varying the radius of the Au inclusion and quantified with the use 
of the filling factor a. 

r a n 

n 16 

Alternatively the volume fraction (f ) , the fraction of conducting to semiconduct- 
ing material, can be defined as follows. 

/=^J=a 2 (2) 

7T7V 

Experimentally disks were grown on a Gallium Arsenide (GaAs) substrate, with 
a 1.3/im thick active layer of InSb. It was noted that the side walls between the InSb 
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Fig. 1. The geometry of the simulated system, illustrating the semiconducting ring (of radius rj,) 
with a conducting inclusion (of radius r a ), including the four point contacts. Here the current is 
flowing through contacts A and B, while the voltage is measured across contacts C and D. 



and the Au regions were approximately 19° from vertical-. The magnetoresistance 
was determined utilising the van der Pauw method of measuring the resistivity 
of a thin filrrP^. This required four Au contacts to be placed equidistant around 
the perimeter of the disk. Two adjacent contacts were used to inject the current 
into the system, with the other two contacts used to measure the resulting poten- 
tial difference. To measure the magnetoresistance the magnetic field was applied 
perpendicular to the surface of the disk. 

In order to simulate this system some aspects were required to be adapted. The 
first simplification made was to simulate the system in two dimensions. This was 
justified since the experimental disk was of a thin film construction. Consequently, 
experimental errors in the angle of the side walls between the two regions were ne- 
glected. Additionally, only the active layer of the disk was considered with all other 
layers of experimental significance ignored for the purposes of simulation. The four 
contacts on the disk were assumed to be point contacts with experimental errors 
disregarded. The simulations were carried out using the finite element analysis soft- 
ware package COMSOL multiphysics. With such assumptions the two-dimensional 
geometry used for simulations is given in Fig. IT] 

Various geometries have been considered, with the filling factor varying by 1/16 
of r;, from to 15/16. With the geometry considered, various parameters were de- 
fined and quantified. The simulations were carried out with an applied field ranging 
from - 5T mirroring those used experimentallySl. In order to measure the mag- 
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netoresistance an electrical current was applied through contacts A and B. Using 
the van der Pauw method for measurement of the resistivity leads to the following 
expression of the magnetoresistance. 

MR _ R(H)-R(0) _ Vcd(H) 

MR - W) ~ ~^W~ " 1 (3) 

Where Vcd(H) = Vc(H) — Vd{H). In order to specify the required mate- 
rials the conductivities and charge carrier mobilities were utilised to establish 
the conductivity tensor for each of the two regions. The conductivity and mo- 
bility values used were as close to the experimental system as possible. The val- 
ues of the conductivity used for InSb and Au were (TinSb = 1-86 x 10 4 (1/$7to) 
and (jau = 4.52 x 10 7 (l/£!m) respectively. With the corresponding mobilities of 
HinSb = 4.55(m 2 /Us) and (j,au =5x 10~ 3 (m 2 /Vs). Using these values, and the 
parameter /3 = fj,H, the conductivity tensor could be defined for each of the two 
materials using the following expression. 

&xx &xy \ _ I J+p2 1+^2 \ ,^, 

Vyx <?yy ) y T+fP 1+/3 2 J 

Finally the boundary conditions were specified. On the perimeter of the entire 
disk (the semiconducting region) they were given as electrical insulation, while at the 
interface between the semiconducting ring and the conducting inhomogeneity they 
were set to continuity. These simulations were carried out with a mesh consisting of 
approximately 60,000 triangular elements, with the mesh being more refined where 
the variation in the potential was greatest. 

3. Results 

Fig. [2] presents the magnetoresistance as a function of applied magnetic field. Here 
we see the magnitude of the magnetoresistance reaching 3,750,000% in a geometry 
with n = 15 (/ = 0.879) and field of 5T. 

The magnetoresistance calculated has an extremely large magnitude already at 
a 5T magnetic field and can be realised at room temperatures. Systems with small 
conducting inclusions (n = 1 - 8) show saturation in the magnetoresistance at a 
magnetic field of approximately 1.5T. Larger inclusions require a larger magnetic 
field to reach their maximum magnetoresistance value and thus saturation occurs at 
higher magnetic fields. The geometries with n = 14 and 15 do not exhibit saturation 
of the magnetoresistance in a field of up to 5T. This indicates that the saturation 
field increases with the size of the conducting region. 

The form and magnitude of the simulated results are in agreement with exper- 
imental data^, with the values especially close at lower applied fields. The magne- 
toresistance being found to be 79% with n = 12 (/ — 0.563) in a field of 0.05T and 
5721% with n = 13 (/ = 0.660) in a 0.25T applied field. 

In a quasi-two-dimensional system it is important that the conducting inclusion 
extends through the entire thickness of the structure. With the application of a 
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Fig. 2. The magnetoresistance as a function of applied magnetic field. Here the symbols corre- 
spond to various geometries: n = 15 (T), 14 (o), 13 (•), 12 (o), 11 ( ), 10 ( ) and 8 (v). The 
largest value of magnetoresistance seen at a field of 5T corresponds to a geometry with n = 15, 
with the lowest value of magnetoresistance at 5T corresponding to a geometry with n = 8. 



magnetic field, perpendicular to the surface of the system (z direction), only charge 
carrier motion in the x-y plane is affected. Charge carriers travel with a motion un- 
altered by the magnetic field when their trajectory is in a parallel direction (along 
the z axis). Therefore, if the conducting inclusion does not extend the entire thick- 
ness of the structure alternate current paths are created (the current can avoid the 
conducting inclusion by travelling underneath). This causes a dramatic reduction 
in the magnitude of the EMR effect. 

The magnitude of the simulated magnetoresistance is found to be larger than 
that found experimentally for n = 14 (/ = 0.766) and 15 (/ = 0.879) especially in 
high fields of approximately 5T. This difference may be a result of the experimental 
system containing side walls between the semiconducting and conducting region 
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Fig. 3. The magnetoresistance as a function of filling factor (a) at 5 various applied magnetic 
fields: H = 0.05T (v), 0.1T (a), 0.25T (▼), IT (a) and 5T (•) 



that deviate from vertical. However, very good agreement is found for smaller Au 
inclusions, for example a magnetoresistance of 39 x 10 4 % for n = 12 (/ = 0.563) 
and 173 x 10 4 % for n = 13 (/ = 0.660) both in a field of 5T. 

Fig. [3] shows the simulated magnetoresistance as a function of filling factor, a. 
This plot shows a good agreement to the experimental results^ yet some differences 
are apparent. We notice that in general the magnetoresistance increases with a 
up to a certain value, with the value of a at which the magnetoresistance peaks 
increasing with higher applied fields. At low magnetic fields (H = 0.05, 0.1 and 
0.25T) the magnetoresistance peaks at n = 12 and for higher fields (H = 1 and 5T) 
peaking at n = 13 and 15. One may notice that the data for the field of 5T does 
not fall after the peak magnetoresistance occurs. Instead the magnetoresistance is 
its largest for n — 15, the geometry with the largest Au inclusion. Experimentally it 
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was found that the magnetoresistance dropped after n — 13 in a 5T field and did not 
continue to increase. This difference may be a result of experimental uncertainties, 
especially those errors in the contacts deviating from ideal point contacts. The 
errors associated with the contacts are more pronounced in geometries where the 
semiconducting ring is narrow (large values of a). 

We represent the path of the total current density for a single system where n = 
12 (the concept generally holds for other geometries) at various applied magnetic 
fields in Figs.[4ta) - (d). The streamlines (red) link points of the same current density 
while the arrows (blue) show the magnitude and direction of the current flow at 
that particular point in the system. These diagrams shows the current injected at 
point A and taken off at point B. 

In Fig.Wa) we see the path of the current flow throughout the system with zero 
applied magnetic field. The current flow here is predominantly directed into the Au 
inclusion, which can be seen most clearly from the arrows (blue), with very little 
of the total current flowing around the perimeter of the disk in the semiconduct- 
ing region. At zero magnetic field the electric field lines are perpendicular to the 
cquipotential Au surface. With no applied magnetic field the current flow is parallel 
to the electric field lines and so flows in the same direction, into the Au inclusion. 
This flow into the conducting region constitutes a system of low resistance. 

In the next diagram Fig.RJb), the applied magnetic field is 0.5T. The current is 
still predominantly flowing through the Au inclusion however the magnetic field has 
had the effect of distorting the current path. Now a larger proportion of the current 
is being forced to flow around the semiconducting region of high resistance. The 
magnetic field causes the electric field (which is still perpendicular to the surface 
of the Au inclusion) and the current flow to stray away from being parallel with 
the electrical field lines. The angle between the electric field and the current (Hall 
angle) increases with higher applied magnetic fields. 

There is a continuation of this trend in Fig. [4?c) where a larger magnetic field 
(H = IT) has caused the angle between the electric field and the current flow 
to increase to near 90°. This means the current flow and electric field are almost 
perpendicular. Thus the current is directed away from the conducting region and 
forced to flow through the semiconducting outer ring of much higher resistance. We 
see here that most of the current (arrow plot) is travelling from point A to point B 
through the outer ring of the system and avoiding the lower resistance Au path. 

Finally, with a high applied magnetic field of H = 5T represented in Fig. E[d) the 
directions of the current flow and the electric field are almost perpendicular. Now 
the majority of the current is flowing around the Au inhomogeneity in the region of 
semiconductor material. The application of the magnetic field causes the current to 
avoid the low resistance path and forces it to flow along narrow (depending on the 
value of a) channels of semiconductor. This results in a system with a substantially 
higher resistance than when no magnetic field is applied. This effect causes the huge 
geometric magnetoresistance known as EMR that was found in experiments^. 

The geometry of the system plays an important role in the magnetoresistance 
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Fig. 4. Visualisations of the path of the total current density in a disk where n = 12 (/ = 0.563). 
The solid streamlines (red) indicate lines of the same current density, while the arrows (blue) 
represent the relative magnitude and direction of the current flow at the given points in the 
system. Figures (a) - (d) show the effect of an applied magnetic field on the flow of current, (a) 
H = 0, (b) H = 0.5, (c) H = 1 and (d) H = 5. 



effect. With a relatively small Au inclusion a large proportion of the current already 
flows through the wide semiconducting channel at H = 0. With the application of 
the magnetic field in such geometries a small proportion of the current switches from 
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(a) Small applied field 



(b) Large applied field 



Fig. 5. Simplified illustration of the current flow in a typical system in both low and high applied 
magnetic fields 



flowing through the Au to through the semiconductor. This leads to smaller values 
of magnetoresistance with low values of a. With large Au inclusions however, more 
of the current initially flows through the Au and thus gets switched to flow through 
the semiconductor when a magnetic field is applied. More of the current is affected 
by the application of the magnetic field, therefore larger values of magnetoresistance 
are observed. Additionally in geometries with large conducting inclusions the semi- 
conducting channels are narrow, this increases the resistance of current flow through 
the semiconductor material, thus further enhancing the magnetoresistance. 

In geometries with large conducting droplets a larger field is required to force 
the current to flow through the narrow semiconducting regions. For this reason 
the geometries with larger values of a peak in their values of magnetoresistance at 
higher applied fields. This explains the trend of the peak magnetoresistance value 
occurring at higher values of a at higher magnetic fields observed both here (Fig. 
pi) and in experimental data 2 . 

These current flow diagrams reinforce ideas discussed by Solin et al. 1 concerning 
the EMR effect. They regard the system as a short circuit (current flowing through 
the conducting region) when the magnetic field is small and an open circuit (current 
deflected around conducting region) when the magnetic field is large. 

The application of the magnetic field can be thought of as a switch; with zero 
or a small applied field the majority of the current flows through the conducting 
Au inclusion giving a low resistance path. However, the application of a large mag- 
netic field forces the current to flow around the perimeter of the disk avoiding the 
conducting inclusion, as if there was a ring of semiconductor with a cavity in the 
centre, thus creating a significantly higher resistance path, see Fig. [5j 
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4. Further Discussion 

The EMR effect strongly depends on the geometry of the system. The physical 
component to the magnetoresistance (ordinary magnetoresistance) being almost 
negligible for both the conducting and semiconducting materials individually. Other 
large magnetoresistance effects originate from a large physical component that relies 
on the electrical properties of the materials involved. The EMR effect may therefore 
be applicable to many different systems, where droplets or regions of material with 
a significantly higher conductivity to that of the surrounding material are observed. 

Due to the high economic value of Au, practical devices based on EMR would 
be more cost effective to produce if the EMR effect was efficient and produced the 
highest value of magnetoresistance with a fixed quantity of Au. In order to optimise 
the amount of Au required the geometry could be modified, for example, with the 
implementation of different shaped droplets (oval shaped droplets instead of circular 
oneeJSl 

Here the effect of a single conducting droplet embedded in a semiconductor is 
investigated, by replicating the experimental geometry in which the EMR effect was 
discovered. However, further research will be carried out investigating the effect of 
multiple conducting droplets embedded in a semiconducting material. Simulation 
of a random droplet model may mean the magnetoresistance follows a linear growth 
with magnetic field as is desirable practically. A multiple droplet system may be 
hard to create precisely but some systems form droplets of conducting and semicon- 
ducting material spontaneously. An example of this is the silver chalcogenides which 
have shown a large room temperature magnetoresistance, where a small excess of 
conductor is added to a semiconducting material. The addition of this small excess 
amount of conducting material has been shown to form a percolating silver net- 
work along the grain boundaries (for already low silver excess) and even conducting 
droplets in the microstructure of the material (for higher excess silver content)^ in 
analogy with this work. It would be interesting to look at a random parquet model 
as a parquet type structure has also been identified^. 

5. Conclusion 

Since this paper has been produced a similar piece of work has come to our atten- 
tion, namely that of Moussa et al\^, which documents finite element simulations 
regarding the EMR effect. Their results show the field dependence of the current 
flow, the potential on the disk periphery and compare the disk resistance with ex- 
periment. These simulations were carried out with a mesh consisting of 6000 nodal 
points. The results given are in excellent agreement with those presented in this pa- 
per and with the experimental data 2 . Our results were calculated using a much more 
sophisticated mesh with the magnetoresistance explicitly calculated and compared 
to the experimental data. 

We have investigated the EMR effect using finite element simulations. The sim- 
ulated geometry was considered in two-dimensions with parameters described in 
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experiments by Solin et al. so as to replicate these systems. The magnetoresistance 
was calculated for different geometries (n = 1 - 15) and in a range of magnetic 
fields (H = - 5T). 

We have found that a huge magnetoresistance appears, due to the geometry of 
the system, and our results show good agreement with the experimental values. The 
size of the conducting inclusion is found to affect the magnetoresistance greatly. We 
have looked at diagrams of the total current density flow throughout the samples at 
various applied magnetic fields and discussed the mechanism for the resulting huge 
magnetoresistance. This model can now be utilised in the investigation of the effect 
of various droplet shapes on the magnetoresistance. 
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